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A.  INTRODUCTION 


The  purpose  of  stand-off  is  survivability  enhancement.  The  price  to  pay 
is  strike  delay.  For  moving  targets,  strike  delay  means  effectiveness  reduc¬ 
tion.  Therefore,  the  effectiveness  analysis  of  stand-off  attack  on  a  moving 
tank  column  must  assess  the  effectiveness  reduction  by  strike  delay.  This  is 
the  prime  subject  of  the  present  methodology  study. 

The  study  addresses  the  following  scenario.  A  tank  column  is  on  the  road 
march  behind  the  FEBA  (Forward  Edge  of  the  Battle  Area).  Most  likely,  but  not 
necessarily,  the  column  moves  toward  the  FEBA.  Blue  has  obtained  reconnais¬ 
sance  about  this  fact  and  plans  a  stand-off  airstrike  or  missile-strike. 

We  use  the  (straight)  road  as  x-axis  and  place  the  FEBA  in  the  positive 
x-direction.  This  has  the  effect  that  the  mean  column  velocity  is  positive. 

The  information  available  to  Blue  may  now  be  described  as  follows: 

1.  At  "reconnaissance  time" 

tr  =>  o  (1-1) 

the  middle  point  of  the  tank  column  occupied  the  position 
x  *  o  (1-2) 

2.  The  average  travel  velocity  of  a  hostile  tank  column  on  the  road  march 
to  the  FEBA  is  <v>.  It  is  assumed  that 


<v>  >  o  (1-3) 

3.  The  length  of  the  column  is  2L. 

4.  The  strike  takes  place  at  the  "strike  time" 

t  >  o  (1-4) 
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The  methodology  offered  in  this  paper  is  designed  to  answer  questions 
of  the  following  type. 

1.  What  is  the  probability  that,  at  strike  time  t,  the  "contemplated 
target  point  T"  with  the  coordinate  x^  is  occupied  by  the  tank  column? 

2.  What  is  the  expected  number  of  tanks  destroyed  if,  at  strike  time 
t,  a  bomb  is  aimed  at  the  aimpoint  A  with  the  coordinate  x_? 

a 

3.  What  is  the  expected  number  of  targets  destroyed  if,  at  strike 
time  t,  m  bombs  are  aimed  at  the  aimpoints  Ap  ...,  Am  with  the 
coordinates  Xp  ...,  xffl  ? 

4.  What  selection  of  aimpoints  maximizes  the  expected  number  of 
tanks  destroyed? 

5.  What  are  the  answers  to  the  aforementioned  questions  if  there  is 
not  one  tank  column,  but  several? 

These  questions  illustrate  only  the  main  thrust  of  the  present  methodology 
study,  but  there  will  be  excursions  and  extensions,  for  example,  a  two- 
dimensional  tank  array  off  the  road,  an  asymmetric  velocity  distribution 
function,  and  others. 

B.  SYMMETRIC  VELOCITY  DISTRIBUTION  FUNCTION 

The  chief  instrument  of  the  analysis  outlined  in  the  previous  section  is 
the  "Velocity  Distribution  Function"  (VDF).  In  the  present  section,  we  assume 
that  the  VDF  is  symmetric  with  respect  to  the  "most  probable  velocity"  7,  also 
called  "mode".  If  nothing  else  is  known  about  the  velocity  distribution,  then 
it  follows  in  conjunction  with  the  Central  Limit  Theorem  of  the  calculus  of 
probability  that  the  VOF  is  Gaussian.  If  ay  denotes  the  "standard 
deviation",  we  have 


f  (v) 


1 

oy  y/2? 


Exp 


iv^L2) 

) 


(2-1) 
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The  Gaussian  VOF  is  constituted  by  two  parameters:  v  and  o  .  These 

parameters  can  be  interpreted  in  terms  of  the  "mean  velocity"  <v>  and  the 

2 

"mean  square  velocity"  <v  >.  To  carry  this  out,  we  consider  the  integral 
equations 

+oo 

r 

f(v)dv  *  1  (2-2a) 

- 


(2-2b) 


+00 

r 


J 

-  00 


f(v)v^dv  3  <v^> 


(2-2c) 


Here,  equation  (2a)  is  simply  the  "normal ization  condition"  which  is  satisfied 
identically.  The  two  remaining  equations  yield 

~  *  <v>  (2-3a) 

*  <v^>  -  <v>^  (2-3b) 

Hence  the  most  probable  and  the  mean  velocity  are  identical.  This,  of  course, 
is  a  consequence  of  the  symmetry  of  the  VOF. 
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Let  us  now  assume  that  a  large  number,  say  N,  observational  or  intelligence 
data  are  available  about  the  motion  of  hostile  tank  columns  on  the  road  march 
to  the  FEBA.  Let  v^  denote  the  travel  velocity  of  the  ith  observation.  We 
then  have 


<v> 


1  n 
■+  5 


vi 


<v1 2>  = 


N 

c 

i=l 


v2 


(2-4a) 


( 2— 4b) 


If  these  equations  are  heeded  in  relations  (3a)  and  (3b),  we  see  that  the 
constitutive  parameters  V  and  can  be  interpreted  in  terms  of  the 
intel 1 igence  data  v. . 

An  example  will  illustrate  this.  We  assume  that  the  following 
intelligence  information  is  available. 

1.  At  any  given  time,  80%  of  the  columns  move  with  approximately  the  same 
speed  of  30  km/h. 

2.  At  any  given  time,  20%  of  the  columns  are  at  rest. 

We  then  have  (with  the  same  approximation): 

Moving:  =  0.8N,  v^  =  30  km/h 

At  Rest:  Nq  =  0.2N,  v^  =  0 
This  yields 

V  =  <v>  =  24  km/h 
cv  =  12  km/h 
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C.  ASYMMETRIC  VELOCITY  DISTRIBUTION  FUNCTION 


Whether  or  not  the  VDF  is  symmetric  depends  on  the  information  which  is 
available  with  regard  to  the  travel  velocities  of  hostile  tank  columns.  We 
shall  elaborate  on  this  subject  in  section  D.  In  the  present  section,  we 
simply  assune  that  the  information  available  dictates  an  asymmetric  VDF.  The 
question  then  is  how  to  construct  this  function. 

Probably  the  most  elegant  way  is  the  following.  We  take  two  Gaussian 
functions 


fx(v) 


,  ,L_ 


Exp 


Iy-v)2 


(3-la) 


f2(v) 


1 


Exp 


(3-lb) 


They  have  the  same  mode,  v,  but  different  standard  deviations  cr ^  anda2> 
For  our  present  purposes,  we  assume 


1  <  a2  (3-2) 

As  in  the  previous  section,  we  denote  the  actual  VDF  by  f(v),  but  it  must 
now  be  kept  in  mind  that  f(v)  is  asymmetric.  The  idea  is  to  "compose"  f(v)  by 
fx(v)  and  f2(v).  For  v<v,  f(v)  shall  be  equal  to  fj(v),  except  for  a 
constant  factor  still  to  be  determined.  For  v>v,  f(v)  shall  be  equal  to 
f2(v),  except  for  another  factor. 

The  purpose  of  the  two  factors  is  two-fold.  First,  at  v*v,  the  two 
functions  must  be  fuzed  together.  And  since 


fx(v) 


._i _ 


f'2(v) 


(3-3a) 
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it  follows  that  equalization  is  achieved  if  f^(v)  is  multiplied  by<^,  and 
f 2( v )  is  multiplied  by  a 

Second,  the  actual  VDF,  f(v),  must  be  normalized.  Therefore,  we  need  one 
free  parameter,  sayay,  to  achieve  normalization.  If  we  then  write  the  two 
factors  in  the  formcj^/j^  and  a2/av,  the  asymmetric  VDF  assumes  the 
following  form 


Figure  1.  Asymmetric  Velocity  Distribution  Function  and  its  Constituents 


The  three  functions  f^(v),  fgCv),  and  f ( v )  are  illustrated  in 
figure  1.  It  can  be  seen  that  f ( v )  is  “biased"  in  favor  of  positive  velocity 
values.  This,  of  course,  is  a  consequence  of  assumption  (3-2)  which,  in  turn, 
accounts  for  our  earlier  assumption  (section  B)  that  the  FEBA  lies  in  positive 
x-di recti  on. 


We  intimated  already  that  f(v)  is  not  yet  normalized,  but  that 
normalization  can  be  achieved  by  proper  selection  of  the  parameter  o  The 
normalization  condition  has  the  form  of  equation  (2-2a).  If  these  equations 
(3)  are  heeded,  we  obtain 


+  00 


i:  f  fi(,)  dv  *  r-_/f2(v> av  ■ 1 


(3-4) 


Here  we  observe  that  f^(v)  and  f ^ (v)  are  normalized  and  that  they  are 
symnetric  with  respect  to  the  axis  v=7.  Hence  it  follows  that 


1 


or 


ay  3  -j-  (jj  +  02)  (3-5) 

Relation  (5)  is  necessary  and  sufficient  for  the  normalization  of  the 
asymmetric  VDF,  f(v).  The  relation  also  demonstrates  that  only  two  of  the 
three  standard  deviations  are  freely  selectable.  If  we 

chose,  for  example,  and  ^  as  freely  selectable,  then  f(v)  is 
constituted  by  three  parameters:  v,  <jj,  02 •  Therefore,  anybody  who 
claims  that  the  VDF  is  asymmetric  should  be  able  to  prove  this  by 
demonstrating  that  the  VDF,  indeed,  is  constituted  by  three  distinct  and 
freely  selectable  parameters.  We  shall  elaborate  on  this  subject  in  the 
following  section. 
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0.  RELATIONS  BETWEEN  THE  CONSTITUTIVE  PARAMETERS  OF  THE  ASYMMETRIC 

VELOCITY  DISTRIBUTION  FUNCTION 

Probabilities  and  distribution  functions  are  based  on  information.  If  the 
information  changes,  so  do  the  probabilities  and  the  distribution  functions. 
The  information  content  of  a  distribution  function  is  stored  in  its 
constitutive  parameters.  However,  in  practical  situations  such  as  stand-off 
attack  on  a  moving  tank  column,  the  information  is  rarely  stored  in  such 
highly  theoretical  paramaters  as  c  ,  Cj,  but  most  1 i kel y  in  more 

observable  parameters  such  as  v,  <v>,  <v2>.  Hence  it  is  desirable,  and  it 
will  also  be  instructive,  to  establish  oy,  o^,  and  as  functions  of  v, 

<v>,  and  <v2>.  This  is  the  purpose  of  the  present  section. 

The  mean  values  <v>  and  <v  >  are  defined  by  equations  (2-2b)  and 
(2-2c).  If  there  we  substitute  relation  (3-3),  we  obtain 


<v>  *  Ij  +  I2 


(4-la) 


<v  >  *  +  J2 


(4-lb) 


V 

•-3J- 


(4-2a) 


U  s  j2-  /  f2(v)vdv 

v  J 


(4-2b) 


V 

*  / 


f^(v)v^dv 


(4-2C) 


J2  *  T~  J  f2(v)y2dV 


(4-2d) 


12 


The  evaluation  of  the  integrals  1^,  I2,  Jp  J2  is  straightforward 
but  tedious.  It  will  be  facilitated  by  the  following  auxiliary  formulae. 
First,  from  equations  (3-la)  and  (3-lb),  it  follows  by  differentiation  that 


vfk(v)  *  V  fk(v)  -  i^- 


(k=l , 2) 


This  may  be  written  in  the  form 


v2fk  *  y2fk  -  ak 


.2  3fk 


(k-1,2) 


For  this,  we  may  also  write 


v2f,,  a  v2f„  -  v  o,2  °^-a} 


k  3 v  uk  37"k 


2~(vfJ  *  l  (k*1. 


2) 


We  also  mobilize  the  following  formulae.  Defining 


E(x)  «  Exp 


x 

7 


and 


■/ 


♦  m  *  /  E(x)  x  dx 


We  have 


<t» 


1 

0  55  T 
\ 


♦l 


2ir 


.  1  ,2 
’2  ~r  A 


(4-3a) 


(4-3b) 


(4-3c) 


(4-4a) 


(4-4b) 


(4-5a) 


(4-5b) 


(4-5c) 
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With  the  aid  of  equations  (4-3a),  (4-3b),  (4-3c),  (4-4a),  (4-4b),  (4-5a), 
(4-5b),  and  (4-5c),  the  integrals  1^,  I2,  Jj,  and  J2  can  now  be 
solved.  The  results  are 


i  ..fL 

1  \  t  2  VzW 


(4-6a) 


02 

o 


_v  +  c2 


•y/77 


(4-6b) 


J,  - 


2a  u 


a  lv 


^  +  3 

—  +al 


yJU' 


( 4-7  a ) 


J2  3 


-2 

J2v  +  — = — 


y]JP 


+ 


(4-7b) 


Now  substituting  equations  (4-6a),  (4-6b),  (4-7a),  and  (4-7b)  into 
equations  (4-la)  and  (4-lb),  we  get 

^(a^+o^) 


<v>  3 


7 


2  2 
o2  -  Oj_ 


(4-8) 


and 


<v2> 


3  3  4v(<72_<71) 

al  +  °2  + 


+  v2  (a^Cj) 


(4-9) 


Under  observance  of  equation  (3-5)  and  the  identity 


3  ,  3 

C1  °2  .  .2  ,2 

°l  •  ala2  +02 


(-T-10) 
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equations  (4-8)  and  (4-9)  assume  the  form 


<v>  =  v 


^fr  (a2-ai> 


<v*>  *  v2  +  2 


£ 


V  (c2  -  Oj)  +  a2-a1a2  +  cr2 


From  equations  (3-5)  and  (4-11),  we  also  obtain  the  useful  relations 


(4-11) 

(4-12) 


For  a  brief  discussion  of  these  results. 
From  relations  (3-2)  and  (4-11)  follows  that 

<  v>  >  *v 


(4-13a) 

(4-1 3b ) 

we  first  observe  the  following. 


(4-14) 


In  other  words,  the  bias  in  favor  of  positive  velocities  places  the  mean 
velocity,  <v>  ,  to  the  right  of  the  most  probable  velocity,!/.  This  result 
follows  also  intuitively  directly  from  figure  1. 


For  further  discussion,  we  solve  equations  (4-11)  and  (4-12)  f or  a ^  and 
This  leads  to  quadratic  equations  which  have  solutions  of  the  form  of 
equations  (13a)  and  13b).  It  is  not  difficult  to  show  that 


v 


<v 


2> 


<  v>  ^  + 


(4-15a) 


Equations  (4-13a),  (4-13b),  and  (4-15a)  establish  oj,  o2,  ov  as 
functions  of  v,  <v>,  and  <v2>.  This  was  the  object  of  the  present  section. 
We  observe  that,  for  "v  ■  <v>,  we  have 


a 

v 


<v> 


T 


(4-15b) 
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Finally,  we  observe  from  relation  (15a)  that  av  is  real  if,  and  only  if. 


-  |  <v>  -  v 

This  condition  is  more  stringent  than  the  condition 

<v2>  -  <v  >2  >  0  (4-16b) 

which  applies  in  the  symmetric  case. 

In  summary,  we  may  now  say  the  following.  The  VDF  is  asymmetric  if,  and 
only  if,  <v>  ^  v.  In  particular,  if  <v>  >  7,  the  VDF  is  biased  in  favor  of 
velocities  greater  than  7.  It  is  frequently  argued  that  we  must  assume  this 
bias  because  of  our  earlier  assumption  that  the  column  moves  in  the  positive 
x-direction.  However,  this  argument  is  faulty  because  the  latter  assumption 
has  already  been  accounted  for  by  the  relation  v  >  o.  Therefore,  the 
assumption  <v>  >  7  constitutes  an  additional  assumption  which  can  only  be 
justified  on  the  basis  of  additional  information.  In  other  words,  the 
symmetry  assumption  <v>  *  7  is  the  analytical  expression  of  the  absence  of 
such  additional  information. 

E.  SYMMETRIC  POSITION  DISTRIBUTION  FUNCTION 

Let  x  denote  the  coordinate  of  the  column's  middle  point  at  strike 
time  t.  This  coordinate  is  not  known  to  Blue  (who  conducts  the  strike),  nor 
is  it  known  to  Gray  (the  analyst).  However,  Blue  and  Gray  have  a  probability 
dP(x,  dx)  that  the  middle  point  lies  in  the  interval  from  x-l/2dx  through 
x+l/2dx.  This  probability  has  the  form 

dP(x,  dx)  =  g(x)dx  (5-1) 

Here  g(x)  is  the  "position  distribution  function"  (PDF).  Since  the  middle 
point  coordinate  is  a  function  of  t,  so  is  the  PDF. 


(4-16a) 
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The  PDF  must  reflect  two  sources  of  uncertainty.  The  first  is  the 
uncertainty  associated  with  the  reconnaissance  information  that,  at 
reconnaissance  time  t  3  o,  the  column's  middle  point  occupied  the  position 
x  3  o.  This  uncertainty  is  customarily  referred  to  as  "Target  Location  Error" 
(TIE).  It  is  also  customary  to  model  the  TLE  by  a  Gaussian  distribution 
function  of  the  form 


g0(xo)  3 


Exp 


(5-2) 


Here,  xQ  is  the  actual  position  of  the  column's  middle  point,  while 


are  the  most  probable  and  the  mean  positions,  all  referring  to  reconnaissance 
time  t  3  o. 

Since  the  distribution  function  go(xQ)  is  "induced"  by  the  TLE 
committed  by  reconnaissance,  we  refer  to  it  as  "Reconnaissance  Induced 
Position  Distribution  Function"  (RPDF). 


The  second  source  of  uncertainty  is  the  uncertainty  about  the  column's 

average  volocity  between  reconnaissance  and  strike.  This  uncertainty  is 

described  by  the  "Velocity  Induced  Position  Distribution  Function"  g(x|xQ). 

This  function  is  subject  to  the  condition  that  the  position  at  reconnaissance 

time  t„  3  o  was  x„. 
r  o 

The  POF,  g(x),  occurring  in  equation  (1)  is  the  convolution  of  the  RPDF, 
90(x0),  and  the  VPOF,  g(xjxQ),  that  is. 


9(x) 


+  00 


/  9(x|Xo>  9o(V  dxo 


(5-4) 
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In  the  present  section,  we  assume  a  symmetric  VDF,  f(v).  Naturally,  the 
VPOF,  gCxjx^),  is  then  also  symmetric.  It  is  now  our  aim  to  derive  this 
function. 

To  this  end,  we  first  establish  the  basic  relations  between  the  column’s 
velocity  and  middle  point  position.  We  have 

x  =  x(t)  =  x0  +  v  t  (5-5) 

Now,  with  respect  to  g(x|xQ),  xQ  is  “given”,  that  is,  constant. 

Therefore,  on  account  of  equation  (5),  the  "velocity  uncertainty"  dv  (not  to 
be  confused  with  the  differential),  induces  the  "position  uncertainty" 

dx  =  t  dv  (5-6) 

It  also  follows  under  observance  of  relations  (3)  that 

x(t)  *  < x ( t )>  =  t<v>  (5-7) 

where  x(t)  and  <x(t)>  are  the  most  probable  and  the  mean  value  of  the  column's 
middle  point  coordinate  at  strike  time  t. 

The  linearity  of  relations  (5),  (6),  and  (7)  implies  that  there  is  a 
biunique  (reversibly  unique)  correspondence  between  the  set  of  all  velocity 
values  and  the  set  of  all  position  coordinates,  given  xQ.  Therefore,  the 
probability  of  a  given  velocity  equals  the  probability  of  the  corresponding 
position.  T'nence  we  have 

dP  |x,  dx  |  x0j  -  dP{  v,  dv}  (5-8) 

Using  here  the  VOF  and  the  VPOF,  this  reads 

g(xjx0)  dx  *  f(v)dv  (5-9) 


10 


If  here  we  substitute  v  and  dv  from  equations  (5)  and  (5),  we  obtain  after 
cancellation  of  dx 


And  if  we  observe  here  equations  (2-1)  and  (7),  we  obtain 


Exp  - 


-  2J 

(x-x  -x(tj) 


2a 


vP 


(5-10) 


(5-11) 


where 


°vP 


to, 


(5-12) 


is  the  standard  deviation  of  the  VPOF. 
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For  a  brief  discussion  of  results  (11)  and  (12),  we  consider  figure  2. 
It  shows  two  trends: 


Figure  2.  Velocity  Induced  Position  Distribution  Functions, 
tj  <  t2,  7  >  o 


1.  With  increasing  strike  time  t,  the  VPDF  travels  to  the  right.  This 
trend  reflects  the  column's  motion  toward  the  FEBA. 

2.  With  increasing  strike  time  t,  the  maximum  of  the  VPDF  decreases, 
while  its  width  increases.  This  trend  reflects  the  increase  of  the  position 
uncertainty  with  increasing  time. 
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Returning  now  to  equation  (4),  we  observe  that  both  g0(xQ)  and 
g(x)xo)  are  Gaussian.  Therefore,  the  convolution  is  also  Gaussian,  anG  the 
result  is 


g(x) 


(jx-x(t)) 


(5-13) 


wnere 


a 


P 


(5-14) 


is  the  standard  deviation  of  the  ?DV. 

For  a  brief  discussion  of  results  (13)  and  (14),  we  observe  that  g(x)  is  a 
function  of  the  strike  time  t  by  virtue  of  the  two  time  functions  x(t)  and 
(t).  Therefore,  g(x)  shows  the  same  trends  as  g(x|xo),  illustrated  in 
figure  Z.  It  also  satisfies  the  normalization  condition 

+<«> 

g(x)  dx  =  1  (5-15) 


F.  ASYMMETRIC  POSITION  DISTRIBUTION  FUNCTION 

If  we  assume  tne  asymmetric  VDF  from  section  D,  the  equations  become  more 
complicated  and,  eventually,  lead  to  two  integrals  which  cannot  be  solved  in 
closed  form  (except  fortfQ  *  o).  However,  the  reasoning  presented  in  the 
previous  section  remains  essentially  unchanged. 

We  first  note  that  equations  (5-1)  through  (5-10)  are  unaffected  by  the 
special  properties  of  the  VDF,  f(v).  However,  equations  (5-11)  and  (5-12) 
must  now  be  replaced  as  follows. 
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g(x|x0) 


5“  9l(xlxo} 


x-x^  x(t) 


x-x  >  x(t) 

o— 


gk(x|x0)  - 


^pV^1 


Exp  /  - 


(x-xQ-x(t))2 

2a2 

<P 


(6-1 


[k  ■  1,2)  (5-2) 


°kp  =  tok  *  1.2) 


(5-3) 


we  then  enter  equation  (5-4)  with  equation  (1),  we  obtain 


a  r 

g(x)  =  ~  /  g1(x|x0)  g0(x0)  dx 

v  x-x(t) 

a 

+  ^  J  92(x|xo)  90(x0)  dxQ 


(6-4) 


The  two  integrals  from  equation  (4)  must  be  solved  either  numerically  or 
approximately.  This  subject  will  not  be  pursued  in  the  present  report.  It 
should  be  noted,  however,  that  for  aQ  *  o,  gQ(x0)  goes  over  into  Dirac's 
Delta  Function.  In  this  case,  of  course,  the  integrals  from  equation  (4)  are 
easily  solved,  and  we  obtain 

Sj (x )  x<x(t) 

g(x)  =  (  (6-5) 

C2 

—  g2(x)  x>  x(t) 

°v 

V 
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with 


gkU) 


(k  ■  1,2) 


(6-6) 


G.  THE  PROBABILITY  OF  OCCUPANCY 

Let  P  (T)  denote  the  probability  that,  at  strike  time  t,  and  before  any 
bombing,  the  "contemplated  target  point"  T  is  occupied  by  the  tank  column. 
The  subscript  zero  stands  for  the  qualifier  “before  any  bombing".  We  shall 
later  also  consider  the  probabilities  of  occupancy  after  the  first  bombing, 
second  bombing,  and  so  forth.  They  will  be  denoted  by  the  symbols  P,(T), 

( T) ,  and  so  forth. 

Figure  3  shows  the  contemplated  target  point  T  with  the  coordinate  xt, 
and  two  possible  middle  point  positions  and  M2  with  the  coordinates 
xt-L  and  xt+L,  where  L  is  the  half-length  of  the  tank  column.  The  points 
and  are  the  lower  and  upper  limits  between  which  the  middle  point 
must  lie  if  the  column  is  to  contribute  to  the  occupancy  of  T.  Therefore, 


Figure  3.  Limits  of  Middle  Point  Positions  Contributing  to  Occupancy  of  T 
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the  probability  of  occupany,  P  (T),  is  the  probability  that  the  column's 
middle  point  lies  in  the  interval  from  x^-L  to  x^+L.  The  latter 
probability  is  obtained  by  integrating  the  PDF,  g(x),  over  the  aforementioned 
interval.  This  yields 


Instead  of  a  copious  discussion  of  the  function  P  (T),  we  offer  figure  4. 
Although  P  (T)  is  not  a  Gaussian  distribution  function,  but  a 
probability,  it  has  the  characteristic  bell-shape.  Its  dependency  on  L  and  t 
is  illustrated  by  considering  two  values  of  each. 


Figure  4.  Probability  of  Occupancy,  v  >  o 
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H.  THE  EXTERIOR  TANK  DENSITY 


I  • 


*> 

{ 


We  assume  that  the  tanks  are  continuously  and  homogeneously*'  distributed 
within  the  column.  If  the  column  has  the  length  2L  and  contains  N  tanks,  the 
number  of  tanks  per  unit  length  is 


If  then  we  denote  the  tank  density  by  p ,  we  have 

In  within 

(3-2) 

o  without 

However,  Blue  and  Gray  do  not  know  the  column's  position  at  strike 
time  t.  Therefore,  they  do  not  know  P  as  a  function  of  x.  But  there  is 
another  tank  density  which  is  known  to  Blue  and  Gray  as  a  function  of  x.  For 
clarity  and  convenience,  we  refer  to  P  as  the  "interior  density",  and  to  the 
density  now  to  be  introduced  as  the  "exterior  density". 

Vie  denote  the  exterior  density  by  the  symbol  v  (x^).  Here  the 
subscript  zero  has  the  same  function  as  at  the  symbol  P  (T),  in  other  words, 
it  is  a  reminder  that  v  (xt)  refers  to  the  state  of  affairs  before  any 
comb i n  g . 

The  exterior  tank  density  is  defined  as  follows.  If  dN  denotes  the  number 
of  tanks  within  the  interval  from  xt-l/2dxt  to  x^+i/2dxt>  then 

dN  »  v  (xfc)  dxt  (3-3) 


1.  It  is  possible  to  assume  a  heterogeneous  tank  distribution,  that  is,  a 
"structured"  column.  'While  this  does  not  pose  any  particular  difficulties,  it 
lengthens  the  report  and  is  therefore  disregarded. 


J 


25 


Since  the  total  number  of  tanks  is  N,  it  follows  that 


w 


*  H 


(3-4) 


Obviously,  the  exterior  tank  density  must  be  related  to  the  probability  of 
occupancy.  Indeed,  if  the  contemplated  target  point  T  is  occupied,  then 
v  (x,.)  =  n.  But  the  probability  that  T  is  occupied  is  P,(T). 

O  W  -yj 

Therefore, 

Vxt}  3  nP0(T)  (3-5) 


If  equation  (5)  is  substituted  into  equation  (4)  and  if,  at  the  same  time, 
equation  (I)  is  heeded,  we  obtain 


/  Pq(T)  dx  t  *  2L  (3-5) 

—  30 


And  if  then  equation  (6-1)  is  observed,  it  follows  that 


/ 


g(x)  dx 


dxt  =  2t 


(3-7) 


We  now  assert 
for  all  values  of 


that  equation  (7)  is  an  identity,  to  wit,  that  it  is  true 
L.  To  prove  the  assertion,  we  introduce  the  special  symbol 


HU  • 


7 


rVL 

Jx,-L 


g(x)  dx 


(8-3) 


The  assertion 


is  proved  if  we  have  demonstrated  that  r(l)  *  21. 
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The  proof  makes  use  of  a  well-known  theorem  of  the  calculus.  Let  the 
function  F(y)  be  defined  as  follows: 

•i>2  (y) 

F(y)  *  f  f (x,y)  dx  (8-9) 

>l(y) 

Here,  f  (x,y),  ip ^  (y),  and  (y)  may  he  any  Function's  subject  only  to 
certain  conditions  of  continuity  and  differentiability.  The  aforementioned 
theorem  then  assumes  the  form 


<o2(y) 

dF  *  f  3f  dx  +£2  f(%  y)  -  lx  ft*  .  y)  (8-10) 

dy  >  (y)  **  dy  ’  dy 


We  apply  this  theorem  to  equation  (8)  for  the  purpose  of  forming  the 
derivative  dT/dL.  Since  equation  (8)  is  a  “repeated  integral",  the  theorem 
must  be  applied  twice.  The  first  application  yields 


dr  * 
dl 


I 


d_ 

dl 


g(x)  dx 


The  second  application  yields 


(8-11) 


dT 

dl 


+  L)  + 


g(xt-  L) 


dx. 


(8-12) 


There  the  substitutions  x^+L  3  5 
and  xt-L  a  Z  yield 


dr  3 

dL 


gU  ) 


•  00 


g(  n)  dn 


(8-13) 
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Now  we  observe  the  normalization  condition  (5-15).  Equation  (13)  then 
assumes  the  form 


d  r  * 

dl 


2 


(8-14) 


It  follows  by  integration  that 

I*  (L)  =  2L  +  const  (8-15) 

If  here  we  set  L  a  o,  we  get 

r  (o)  a  const  (8-16) 

Hence  we  have 

r ( L) *  2L  +  r (o)  (3-17) 

But  from  equation  (8)  can  be  seen  that  r (o)=o.  Hence,  it  follows  that 

r  (l)  a  21  (8-18) 

which  was  to  be  proved. 
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I.  ATTACK  GEOMETRY 


Figure  5  illustrates  the  attack  geometry.  It  shows  the  two  Cartesian 
coordinate  systems  (x,  y)  and  (E,n)  with  all  axes  parallel  to  the  ground.  The 
three  points  A,  I,  T  and  their  coordinates  are  defined  in  table  1. 


Table  1 


DEFINITION  OF  POINTS  AND  COORDINATES 


NAME 

SYMBOL 

COORDINATES 

Aimpoint 

A 

V  V  V  na 

ra 

Impact  Point 

I 

V  yi*  ?i’  ni 

7. 

i 

Target  Point 

T 

xt’  yt*  5t"  nt 

7t 

As  before,  the  road  on  which  the  tank  column  moves  toward  the  FEBA  serves 
as  x-axis.  The  y-axis  serves  to  measure  the  offset  from  the  road.  The 
trajectory  of  the  attack  vehicle  (aircraft  or  missile),  serves  as  E-axis. 
Then-axis  measures  the  offset  from  the  trajectory  of  the  attacker.  The  angle 
between  x-axis  and  E-axis  is  denoted  by  3  and  called  the  "attack  angle". 

The  contemplated  target  point,  T,  may  lie  anywhere  on  the  road.  The 
reason  is  that,  with  every  point  on  the  road  (except  the  infinities),  a  non- 
vanishing  probability  is  associated  that  it  may  be  occupied  by  the  colunn. 
Therefore,  every  point  on  the  road  is  a  "potential"  target  point.  However, 
since  the  tanks  are  confined  to  the  road,  it  follows  that 

yt  -  o  (9-1) 

Since  the  potential  target  points,  T,  are  confined  to  the  road,  so  are  the 
aimpoints.  A,  that  is, 

ya  *  o  (9-2) 

The  impact  point,  I,  is  randomly  related  to  the  aimpoint,  A.  Thence,  it 
has  a  generally  non-vanishing  offset,  y1 ,  from  the  road. 

The  radius  vectors  ra,  r^,  and  r*,  shown  in  table  1,  but  not  in 
figure  5,  refer  to  the  origin,  0,  of  the  (x,  y)  coordinate  system.  Shown  in 
figure  5  are  the  "distance  vectors" 


rti  3  ri  *  rt 


rta  3  ra  '  rt 


r.  3  r  -  r. 
la  a  l 


(9-3a) 

(9-3b) 

(9-3c) 
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They  point,  in  this  order,  from  T  to  I,  from  T  to  A,  and  from  I  to  A.  For 
the  norms  of  these  vectors,  we  obtain  under  observance  of  equations  (1)  and 
(2)  the  relations 


2  a  (x  .  x  )2  +  y2 

ti  'xi  xt‘  yi 

(9-4a) 

,2  ♦  (x  -  x,)2 

ta  5 

( 9— 4b ) 

.2  -  («,  -  «1)2  *  y(2 

i  a 

(9 -4c) 

From  figure  5,  we  read  the  following  equations  of  coordinate 
transformation: 


x  3  x.  +  5  cos  8  -  n  sin  S 

a 

(9-5a) 

y  3  5  sin  6  +  n  cos  B 

(9-Sb) 

5  3  (x-xj  cos  8  +  y  sin  8 

a 

(9-6a) 

n  *  -( x-x  )  sin  8  +  y  cos  8 

a 

(9-6b) 

THE  AIMPOINT  RELATED  KILL  PROBABILITY 

The  ordinary  kill  probability  is  "impact  point  related". 

It  has  the  form 

Pk<T|l>  ■'(■•«) 

(10-1) 

This  is  the  probability  that  a  target  at  T  is  killed  by  a  munition 
impacting  at  I.  Generally,  this  probability  is  a  function  of  the  distance 
rfci  from  T  to  I,  as  Indicated.  In  the  present  report,  we  consider  the 
impact  point  related  kill  probability,  that  is,  the  function  f(rti),  as 
given  by  sources  external  to  this  report. 
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The  object  of  the  present  report  is  the  "aimpoint  related  kill 
probability"  Pk(T|A).  This  is  the  probability  that  a  target  at  T  is  killed 
by  a  munition  aimed  at  A.  For  the  calculation  of  this  probability,  we  neeo 
the  probability  that  a  munition  aimed  at  A  will  impact  at  I.  Denoting  this 
probability  by  dP ( I J A) ,  we  have 

Pk(T|A)  3  j*p’K(T\V  dP(l|A)  (10-2) 

I 

Here,  the  integration  is  to  be  extended  over  all  possible  impact  points,  I. 

As  usual,  we  generate  the  probability  dP(ljA)  by  an  "Impact  Distribution 
Function"  (IDF).  Denoting  the  IDF  by  the  symbol  h( I | A) ,  we  have 

dP ( I  (A)  =  h ( I JA)  dXi  dy.  (10-3) 

Substituting  equation  (3)  into  equation  (2),  we  obtain 


+oo 


-00 


(10-4) 


If  here  the  functions  pk(T|I)  and  h( 1 1 A)  a r/»  known,  the  aimpoint  related 
kill  probability  Pj^ ( T | A)  can  be  calculated  on  the  basis  of  this  equation. 


Strictly  speaking,  the  IDF  lies  outside  of  the  scope  of  this  report. 
However,  since  it  is  convenient  and  simple,  we  may  as  well  provide  this 
function  at  this  occasion.  We  first  note  from  figure  5  that  the  trajectory  of 
the  attacker  has  S-direction.  Hence,  we  use  for  the  formulation  of  the  IDF 
the  (  5,  n)-coordinate  system,  switching  to  the  (x,y)-coordinate  system 
thereafter.  Also  noting  from  figure  5  that  the  aimpoint.  A,  has  the 
coordinates 
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Here,  P  (T)  and  P,  (T|a  )  are  given  by  equations  (7-1)  and  (10-4),  The 

0  K  I 

corresponding  exterior  tank  density  is 


v1  (xt)  *  nP]  (T) 


(11-2) 


This  is  the  tank  density  at  the  arbitrary  point  T  after  one  bomb  has  been 
aimed  at  Ap 

After  m  bombs  have  been  aimed  at  Ap...,  Am,  the  probability  of 
occupancy  of  T  is 


p»  (T>  *po  (T)  TT  i’  -  "k  '! 

a  «  1  ' 


(11-3) 


and  the  corresponding  exterior  tank  density  is 


Um  {xt)  =  nP_(T) 


(11-4) 


Let  then  <  ;,,s JAi *  •••»  Am>  and  <  I A1 . A'->'  rsspect1vely, 

denote  the  expected  numbers  of  tanks  that  have  survived  and  that 

have  been  killed  after  m  bombs  have  been  aimed  at  Ap  ...»  A^.  If  it  is 

then  recalled  that  the  total  number  of  tanks  in  the  column  is  N,  it  follows 


<n$K . v  * <  \lAi . v  • N 

On  the  other  hand,  we  also  have 


(11-5) 


<NS|A1  .  in’  '  X 


Hence,  it  follows  that 


(11-6) 


<  Nk  I  Ap  ....  V  *  N  -  n  J  (T)  dxt 

I  •  00 


(11-7) 


If  here  eauation  (3)  is  observed,  we  obtain 


<  ^k  A1  *  *  *  * » 


/  Po(T’  aVi  }  1_Pk  (T  |  Aa)  j  dxt  {U-8) 


*  N  -  n 


r 

’a 


=  o 


n.  *  o 

a 

we  write 


h(Itfl)  -55^-  ExP 


(10-Sa) 

(10-Sb) 


(10-6) 


Here  and  are  the  standard  deviations  in  ''range"  and 
"deflection",  that  is,  in  5-direction  and  in  n-direction.  The  corresponding 
"errors  probable"  are  related  to  the  standard  deviations  by  the  equations 


REP  =  0. 6740 ^ 


(10-7a) 


DEP  =  0.674°0 


(10-7b) 


For  the  transformation  from  the  (  5,  n)-coordinates  to  the 
(x, y)-coordinates,  we  infer  from  equations  (9-6a)  and  (9-6b)  that 


?i  *  (xi**a)  cos  6  +  yi  sin  6 

(10-8a) 

n_-  *  sin  9  +  y.  cos  9 

lid  1 

(10-8b) 

If  these  equations  are  substituted  into  equation  (6),  the  IDF  ■ 
the  form  in  which  it  is  needed  in  equation  (4). 

is  given  in 

K.  EXPECTED  NUMBER  OF  TANKS  KILLED 

We  recall  from  section  G  that  P  (T)  denotes  the  probability  of 
occupation  of  the  contemplated  target  point  T  before  any  bombing. 

P 1 (T)  denote  the  probability  of  occupation  of  T  after  one  bomb  has 
aimed  at  A.  We  then  have 

Let  now 

been 

P1  (T)  •  P0(T)  j  1  -  Pk(T|A  )  j 
{  *34 

(11-1) 

(11-1) 


For  computational  purposes,  it  is  convenient  to  conver:  the  product  under 
the  integral  into  sums  of  products,  applying  the  formula 


rr 

(1-x.)  = 

m 

i  -E 

m-1 

x  +  C 

i  =  1 

l 

i  *  1 

i»l 

m-2 

m-1 

m 

c 

E 

E 

x.  xk  ^  +  - 

i=l 

k  =  i+l 

l  sk+l 

.X.  - 
1  k 


(11-9) 


When  applying  this  formula  to  equation  (8),  we  observe  in  conjunction  with 
equations  (8-1)  and  (8-6)  that 


+® 


N  -  "  /  P„(T)  dx,  * 


(11-10) 


Equation  (8)  then  assumes  the  form 

<Nk|Al . Am>  =n  /  PQ  (t)  C  Pk  (T  |  Ai)  dx 

•oo  i=l 

•"  2  po(t>  H  4,i  p<  <t  iai>  dxt*  ■  *  ••• 


(li-ii) 


Equation  (11)  is  the  general  result  for  arbitrarily  selected  aimpoints 
Al»  ...»  Am<  For  a  brief  discussion  of  this  result,  we  consider  the  case 
m  =  1.  Equation  (11)  then  reduces  to 


+<» 

Nk |  A >  -  n  J  p0(T)pk(T|A)dxt 


(11-12) 
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If  here  equations  (7-1),  (8-1),  and  (10-4)  are  heeded,  this  becomes 


400 


<  N,  A>  = 


VL 


J 


g(x)dx 


VL 


4.  co 

j*  Pk( Tj I)h(l|A)dx.  dy. 


dx„ 


(11-13) 


Here  we  offer  the  following  observations: 


1.  The  aimpoint,  A,  is  freely  selectable.  However,  it  is  not  difficult 
to  see  that,  for  a  symmetric  PDF,  g(x),  the  optimal  aimpoint  is  the  most 
probable  position  of  the  column's  middle  point,  that  is, 

xa  =  x(t)  (11-14) 

Since  >T(t)  =  t7,  it  follows  that  the  optimal  aimpoint  travels  with  the 
column's  most  probable  velocity  toward  the  FEBA. 

The  optimal  aimpoint  for  the  asymmetric  PDF  will  be  determined  in  section  L. 

2.  By  virtue  of  the  PDF,  g(x),  the  expectation  value  <NjA>  depends  on 
the  strike  time,  t.  More  specifically,  <N ^ j A>  decreases  as  t  increases. 

3.  Also  by  virtue  of  g(x),  <Nk|A>  depends  on  aQ  and  ay,  the  standard 

deviations  associated  with  the  TLE  and  the  VDF.  More  specifically,  <Nk|A> 

decreases  as  a  and  a  increase, 
o  v 

4.  From  equation  (13)  can  be  seen  that  <Nk|A>  depends  on  the  column's 
half-length,  L.  It  can  be  shown  that  <N|JA>  decreases  as  L  increases, 
provided  N  is  kept  constant. 

5.  By  virtue  of  the  impact  point  related  kill  probability,  Pk (T( I ) , 

<Nk|A>  depends  on  the  kill-potential  and  the  kill-range  of  the  munition. 
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6.  By  virtue  of  the  IDF,  h(l|A),  <NjA>  depends: 

a.  On  the  aimpoint  selection. 

b.  On  the  attack  angle,  e. 

c.  On  the  delivery  accuracy  measured  by  oR  and  <NJA> 
decreases  as  aR  and  aD  increase. 

L.  OPTIMAL  AIMPOINT  FOR  ASYMMETRIC  POSITION  DISTRIBUTION  FUNCTION 

We  consider  the  case  m  =  1,  that  is,  one  bomb  per  tank  column.  In  case 
of  the  symmetric  PDF,  it  is  obvious  that  the  optimal  aimpoint  is  identical 
to  the  most  probable  middle  point  position,  iT(t),  which  in  turn,  is  identical 
to  the  mean  middle  point  position,  <x"(t)>.  However,  in  case  of  the  asymmetric 
PDF  (3-3),  the  mean  middle  point  position  lies  somewhat  to  the  right  of  the 
most  probable  position.  For  similar  reasons,  it  is  to  be  expected  that  the 
optimal  aimpoint,  too,  lies  somewhat  to  the  right  of  >f(t).  This  is  illustrated 
in  figure  6,  where  A*  denotes  the  optimal  aimpoint  with  the  coordinate  x*  . 

a 


Obviously,  the  condition  for  aimpoint  optimization  is 
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This  equation  must  have  a  unique  solution  for  which  we  denote  by 
xa<  We  then  recall  from  section  6  that  P  (xt)  is  generated  by  integra¬ 
tion  of  g(x)  over  the  interval  from  xt  -  L  to  x{  +  L.  Since  g(x)  is  now 
asymmetric  with  respect  to  x(t),  the  question  arises  how  3T(t)  relates  to 
the  interval  of  integration.  In  figure  6,  we  have  tacitly  assumed  that 
x(t)  lies  within  the  interval,  that  is. 

Case  I:  xt  -  L  £  x(t)  <  x t  +  L  (12-2) 

Principally,  however,  we  also  have 
Case  II:  x"(t)  <  x^  -  L 
Case  III:  >T(t)  >_  x^  +  L 

One  can  show,  however,  that  Cases  II  and  III  do  not  yield  a  maximum  of 
PQ(xt ) -  This  can  also  be  inferred  intuitively  from  figure  6.  In  order  to 
save  time  and  space,  we  proceed  immediately  to  Case  I. 

If  condition  (2)  is  satisfied,  it  follows  from  equations  (7-1),  (3-la), 
and  (3-lb)  that 


xt  +  L 

g j(x)  dx  ♦  °Z  f  g2(x)dx  (12-3) 
av  J-(t) 


x(t) 


u 
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To  form  the  derivative  3PQ  |  3xt,  we  utilize  the  auxiliary 
theorem  from  section  H.  This  yields 


g1(xt-L)  J-  K-L>  -  ll  g^xt+i-)  3_  (x.h)  (12-4) 


If  here  we  set  3PQ  /  3x^.  =  0,  we  must  replace  xt  by  x*a.  This 
yields 


al  9i  (x*  -L)  =  f_2  gJx*  *  L) 
a  1  a  a  L  a 

v  v 


(12-5) 


Here  we  observe  equation  (6-6)  and  obtain 


(xa  -  L  -  x  ( t) ) 

7^- 


(xa  +  L  -  x  ( t ) ) 

V 


(12-6) 


Taking  the  square  root,  we  must  make  sure  that  both  expressions  are 
positive  in  accordance  with  condition  (2).  This  yields 


x*a  -  L  -  x  (t)  ,  x  (t)  -  x*a  -  l 

7[  ~2 

This  equation  is  linear  in  x£  and  has  the  solution 

x\  ■  xU)  +  L  a2  -  CT1 
3  *0, 


(12-7) 


(12-8) 
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This  is  the  optimal  aimpoint  coordinate. 

For  a  brief  discussion  of  this  result,  we  first  consider  the  symmetric 

★  — 

case  =  aj.  It  yields  xa  *  x(t),  as  expected.  Next  we  consider 

the  case  L  s  o.  It  also  yields  x  »  x(t),  as  should  be  expected. 

®  *  — 

Finally,  we  recall  condition  (2).  It  yields  x.  >  x(t),  as  expected. 

cl  ' 

M.  DISTRIBUTION  FUNCTIONS  FOR  RELATIVE  VELOCITY  ANO  SPACING  BETWEEN  COLUMNS 

In  the  next  section,  we  shall  consider  an  aggregate  of  more  than  one 
column.  The  prerequisites  for  this  problem  are  the  “Relative  Velocity  Distri¬ 
bution  Function11  (RVOF),  and  the  “Spacing  Distribution  Function"  (SDF),  which 
will  now  be  der i ved . 


Let  dp  |  (v,  dv);  (v-vr,  dvf)  j  denote  the  joint  probability  that  one 
column  has  a  velocity  in  the  interval  from  v  -  l/2dv  to  v  +  l/2dv,  and  the 
other  column  has  a  velocity  in  the  interval  from  v-v  -  1/2  dvr  to  v-vr 
+  l/2dvf.  We  then  have 


dP  j  (v,  dv);  (v-vr,  dvr)  j  *  f(v)  dvf(v-vr)  dvr  (13-1) 

Equation  (1)  describes  one  case  of  an  infinite  set  of  cases  in  all  of 

which  the  two  columns  have  the  relative  velocity  v  Let  dP(vr,  dvr) 

denote  the  probability  of  the  total  set,  that  is,  the  probability  that  the 

relative  velocity  between  the  two  columns  lies  in  the  interval  from  v„  - 

r 

l/2dv  to  v  +  l/2dv„.  Introducing  the  RVDF,  f  (v  1,  we  write 

r  r  r  p  P 


dP(vr,  dvr)  *  f  (vr)  dvr 


(12-2) 
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From  the  definition  of  dP(vr,  dvr)  follows  that  it  is  obtained  by 
integration  of  equation  (1)  over  all  values  of  v  from  -»  to  +08  .  This  yields 


dP(vr,  dvr) 


f(v)  f(v-vr)  dv 


(13-3) 


Observing  here  equation  (2),  we  obtain 

+00 

f  (v  )  =  f  f( v )  f(v-v  )  dv 

r  r  J 


(13-4) 


In  other  words,  the  RVDF  is  the  convolution  of  the  VDF  with  itself.  We 
may  also  say  that  it  is  the  "autocorrelation  function". 

If  f(v)  is  Gaussian,  so  is  fr(vr).  Assuming  now  the  symmetric  VDF 
from  section  B,  we  have 


f  ( v  ) 
r  r 


(13-5) 


with 


%  *  %  \[7^  U3-6) 

For  a  brief  discussion  of  this  result,  we  first  observe  that  the  most 
probable  and  the  mean  value  of  the  relative  velocity  are 
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Vr  *  <vr>  =  0 


(13-7) 


as  should  be  expected.  Next  we  observe  that  a  >  a  .  This  relation 

r  v 

reflects  the  fact  that  two  independent  col  wins  constitute  a  system  of  two  sta 
tistical  degrees  of  freedom.  Naturally,  such  a  system  leaves  more  room  for 
error  than  a  system  with  only  one  degree  of  freedom. 

Next  we  turn  to  the  "Spacing  Distribution  Function”  (SDF).  Let  sQ  and  s 
denote  one  half  of  the  distances  between  the  middle  points  of  two  adjacent 
columns  (half-spacing)  at  reconnaissance  time  t  *  o  and  at  strike  time  t. 

We  then  have 


s  »  sQ  +  vrt  (13-8) 

and 

ds  3  tdvr  (13-9) 

Let  then  dP(s,ds)  denote  the  probability  that  the  spacing  lies  between  s 
-  1/2  ds  and  s  +  1/2  ds.  Denoting  the  SDF  by  F(s),  we  write 

dP(s,  ds)  3  F(s)  ds  (13-10) 

We  then  infer  from  the  linearity  of  equations  (3)  and  (9)  that 

dP(s,  ds)  *  dP(vr,  dvp)  (13-11) 
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Hence  we  have 


F(s)  ds  =  fr(vr)  dvr  (13-12) 

Here,  we  observe  equations  (8)  and  (9).  This  yields 

F(s)  ■  1  fr  (12-13) 

And  this  yields 

F(S)  *  _ 1_ 

as 

with 

*  tar  (13-15a) 

or 

%  -  yfT  to  v  (13-lSb) 

For  a  brief  discussion  of  this  result,  we  first  observe  that  the  most 
probable  and  the  mean  value  of  s  are  identical  to  the  original  reconnaissance 
estimate  so,  that  is. 


Exp  l  -  (_s-so) 


(13-14) 


2  a. 


s  s  <s>  *  s  (13-16) 

This  is  in  agreement  with  relation  (7).  Next,  we  observe  that  a$ 
increases  with  time.  This  simply  reflects  the  general  rule  that  a  velocity 
error  leads  to  a  time-proportional  position  error. 
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N.  AGGREGATES  OF  MORE  THAN  ONE  COLUMN 


If  there  are  more  than  one  column  on  the  road,  the  probability  of  occu¬ 
pancy,  P  (T),  is  increased  because  each  column  is  a  potential  contributor 
to  the  occupancy  of  T.  We  now  assume  that  reconnaissance  has  detected  an 
"aggregate"  of  M  columns  on  the  road  march  toward  the  FEBA.  For  convenience, 
we  assume  that  M  is  odd.  In  this  case,  the  aggregates's  middle  point  is  the 
middle  point  of  the  middle  column.  At  reconnaissance  time 
t  =  o,  this  middle  point  occupies  the  position  x  »  o. 

Also  for  convenience,  we  assume  that  all  columns  are  equal.  This  implies 
that  they  have  the  same  half-length,  L,  that  they  contain  the  same  number  of 
tanks,  N,  and  -  most  important  -  that  they  obey  the  same  VDF,  f(v). 

The  challenge  of  the  present  problem  lies  in  the  task  of  providing  a 
realistic  description  of  the  interdependency  of  the  columns  in  analytical 
terms.  Certainly,  the  columns  may  have  statistically  different  velocities, 
though  they  obey  the  same  VOF,  but  they  must  not  penetrate  each  other  or 
change  their  order.  Hence,  they  are  not  quite  independent  of  each  other,  but 
they  are  not  rigidly  connected,  eithef 

For  the  analytical  treatment  of  this  problem,  we  derived  in  the  previous 
section  the  SDF,  F(s).  From  the  derivation  of  this  function  follows  that  it 
assumes  complete  independence  of  the  columns.  This  can  also  be  inferred  from 
the  fact  that  s  is  allowed  to  assume  any  value  including  negative  values. 

This,  however,  violates  the  condition  that  the  columns  must  not  penetrate  each 
other  or  change  their  order.  This  condition  will  now  be  enforced  by 
"renormalization"  of  the  SDF. 

We  recall  that  L  denotes  the  half-length  of  the  single  column,  and  that  s 
denotes  the  half-spacing,  that  is,  half  of  the  distance  between  column  middle 
points.  Using  these  two  concepts,  we  now  express  the  condition  of  no  pene¬ 
tration  by  the  following  inequality: 


s  >  L 


(14-1) 


U 


This  condition  is  satisfied  if  we  replace  the  SDF,  F(s),  by  the  following 
function 


R(t)  F(s)  s>l 

F„(s)  -  (14-2) 

K  o  s  <  L 

Here,  R{t)  is  the  "renormalization  function"  which,  as  indicated,  depends 
on  the  time.  It  is  the  solution  of  the  "renormalization  condition" 


R(t)  /  F(s)  ds  =  1 
L 


(14-3a) 


which,  in  view  of  equation  (2),  may  also  be  written  in  the  form 


FR(s)ds  =  1  (14-3b) 

The  function  FR(s)  is  called  "Renormalized  Spacing  Distribution  Function" 
(RSDF).  It  satisfies  the  condition  of  no  penetration  while  imposing  no 
further  conditions  on  the  free  mobility  of  the  column. 
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Figure  7.  Aggregate  of  Five  Columns 
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Figure  7a.  Arrival  and  Departure  of  First  Column 


Figure  7b.  Arrival  and  Departure  of  Second  Column 


Figures  7a  and  7b  show  an  aggregate  of  five  columns,  numbered  1  through  5 
in  the  order  in  which  they  pass  through  any  given  point.  The  aggregate's 
middle  point  is  denoted  by  M.  In  the  present  example,  this  is  the  middle 
point  of  the  third  column. 

In  both  figures,  the  spacings  between  the  middle  points  of  adjacent 
columns  are  denoted  by  2sj ,  2s2,  2S3,  and  2s^.  These  are  independent  random 
variables  obeying  the  same  RSDF. 

In  both  figures,  attention  is  focused  on  the  contemplated  target  point  T 
with  the  coordinate  x^.  The  upper  and  lower  bars  connect  the  five  columns 
in  different  "situations".  In  figure  7a,  the  upper  bar  connects  the  columns 
in  a  situation  in  which  the  first  column  just  "arrives"  at  T,  while  the  lower 
bar  depicts  the  situation  where  the  first  column  "departs"  from  T.  In  figure 
7b,  the  "arrival"  and  "departure"  depicted  by  the  upper  and  lower  bars  refer  to 
the  second  column. 

It  is  important  to  understand  that  the  four  situations  depicted  in  figures 
7a  and  7b  do  not  necessarily  refer  to  four  different  times.  At  any  given  time 
each  of  the  four  situations  is  possible.  Of  course,  the  situations  are  mutu¬ 
ally  exclusive,  that  is,  only  one  of  them  is  possible  at  any  given  time. 

To  each  of  the  four  situations  corresponds  a  position  of  the  middle  point 
of  the  aggregate.  In  figure  7a,  the  positions  corresponding  to  the  upper  and 
lower  barsare  x^  and  x^.  and  in  figure  7b,  they  are  x21  and  x22.  Hence  the 
first  index  indicates  the  column,  while  the  second  index  indicates  arrival  or 
departure. 

We  then  read  from  figure  7a  the  following  relations: 

x2 1 ( S1 ,  s2)  +  2(sj  +  s2)  -  L  =  xt  (14-4a) 

X1 2 ( S1  ’  S2^  +  2^sl  +  s2^  +  L  “  xt  (14— 4b) 
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Here  the  notation  x^U^,  s2)  and  x12(Sj,  s2)  indicates  that  these  positions 
are  functions  of  s^  and  s2.  For  convenience,  we  combine  the  two  equations  (4a) 
and  (4b)  in  the  form 

xlk^sl’  s2^  +  2^sl  +  s2^  -  L  =  xt^k  =  lf  2^  (14-5a) 

where  the  upper  and  lower  signs  of  L  refer  to  the  upper  and  lower  bars  of 
figure  7a  or,  what  amounts  to  the  same,  to  k  =  1  and  k  +  2,  indicating  arrival 
and  departure. 

In  a  similar  way,  we  obtain  from  figure  7b  ' 

x2k(s2)  +  2s2  +  L  -  xt(k  =  1,  2)  (14-5b) 

where  the  first  index  of  x2k  indicates  that  now  the  second  column  is  under 
consideration.  If  then  we  consider  also  the  third,  fourth,  and  fifth  columns, 
we  obtain 

x3k  ±  L  =  xt(k  =  1 »  2)  (14-5c) 

x4k(s3)  -  2s3  +  L  -  X t ( k  =  1,  2)  (14-5d) 

x5k(s3,  s4)  -  2(s3  +  s4)  +  L  =  xt(k  =  1,  2)  (14-5e) 

The  middle  point  positions  xlk(s1,  s2),  s2k(s2),  x3k,  x4k(s3),  and  x5k(s3,  s4) 
are  subject  to  the  condition  that  the  spacings  assume  the  values  listed  in 
parenthesis.  Therefore,  if  we  integrate  the  PDF,  g(x),  over  the  corresponding 
intervals,  for  example,  from  x^(s^,  s2)  to  -y.j2(Sj,  s2)»  we  obtain  the  corres¬ 
ponding  conditional  probabilities  of  occupancy.  Me  write 

xl2(si>  s2) 

P1(xt|si,s2)  = 

X11 (S1  ’  s2* 


/ 


g(x)dx 


(14-6a) 
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*22^ s2^ 


P2(x 


tU2)-  / 


g(x )dx 


x£i ( s2) 


'32 


p3(xt)  =  g(x)dx 


31 


x42(s3) 


P4(xt|s3>  =  /  g(x)dx 


x4l(s3) 


(14-6b) 


(14-6c) 


( 1 4  -  6d ) 


x 52 ^ s3 *  S4^ 
P5(xtlS2*  s4>  =  j  g(x)dX 

x 51 ( S3 ’  S4^ 


(14-6e) 


Me  then  proceed  from  the  conditional  to  the  unconditional  probabilities 
of  occupancy  by  convolution  with  the  RSDF.  This  yields 


+» 


Pl(xt>  -JT  P1(X,|S1-  s2>  FRlsl>  FR(s2>dsl  ds2 


(14-7a) 


+00 


P2(xt) 


J  P2(xt|s2)  FR(s2)ds2 


(14-7b) 


P3(xt)  “  P3(xt} 


(14-7c) 
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+CO 


P4(xt}  =/P4(xt|s3}  FR(s3)ds3 


(14-7d) 


■foo 


P  5  ( x  t )  = 


j j  P5^Xt^4’  7  FR^s4^  FR^s5^ds4  ds5 


(14-7e) 


We  now  make  certain  simplifications  which  are  significant  in  terms  of 
the  reduction  of  volume  and  labor,  but  rather  insignificant  in  terms  of  the 
errors  which  they  introduce.  To  that  end,  we  observe  from  equations  (5a) 
through  (5e)  that  the  s.  occur  only  as  sums,  for  example,  s^  +  Sg.  Therefore, 
equation  (7a),  for  example,  has  the  form 


+oo 


w 


ff  *(S1  +  S2}  FR(s1)  FR(s2)dsl  ds2 


(14-8) 


We  assert  that  this  double  integral  can  be  approximated  by  the  following 
single  integral: 


+oo 


pl(xt 


)  = 


7 


♦  (2s) 


FR(s)ds 


(14-9) 


Since  the  proof  of  the  assertion  is  rather  involved,  it  has  been  relegated 
to  the  appendix. 
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If  these  simplifications  are  carried  out  in  equations  (5a)  through  (5e), 
we  obtain 


Xlk  (S)  s  xt  -  4S  l  L  (14-10a) 
X2k  (S)  s  xt  *  2S  +  L  ( 14— 10b ) 
X3k  «  Xt  l  L  f  (  k  =  1,  2)  (f4-10c) 

X4k ( s)  *  Xt  +  2S  ;  L  (14-I0d) 
X5k  (S)  =  Xt  +  4S  ;  L  J  (14-I0e) 


Equations  (6a)  through  (6e)  then  assume  the  following  form 


X]^(  s) 

pi  ^Xt  j  =  /  g(*)  dx 

X11 ( S) 


P2  (Xt  S) 


X22(S) 

3  J  g(x)  dx 
x71(S) 


X12(S) 

P3  (Xt  |  s)  *  f  g(x) 
Xj1(S) 


P4  (Xt  S)  * 


X,,(S) 


x4i(s) 


g  (x)  dx 


( 14-11  a ) 


( 14-llb ) 


( 1 4—11 c ) 


( 1 4-1 Id ) 
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p. 

3 


(X, 


S)  * 


x„(S) 

7 

*51  ‘SI 


9  (X)  dx 


(14-lle) 


Equations  (6a)  through  (6e)  establish  the  conditioned  probabilities  that 
the  first  through  fifth  columns  occupy  T.  3ecause  of  the  condition  of  no 
penetration,  these  prooabil ities  are  mutually  exclusive.  Therefore,  the  prob¬ 
ability  that  any  one  of  the  five  columns  occupies  T  is 


5 


E 

i=l 


(14-12) 


We  then  proceed  from  the  conditional 
occupancy  by  convolution  with  the  RSDF. 


to  the  unconditional  probability  of 
This  yields 


<V 


+oo 

/ 


P0  (Xt  S)  Fr  (S)  ds 


(14-13) 


Equations  (10),  (11),  (12),  and  (13)  permit  the  calculation  of  the  joint 
probability  of  occupancy  of  any  point  T  on  the  road  at  any  time  by  the  columns 
of  an  aggregate  of  five.  Extension  to  aggregates  of  M  columns  is  straight¬ 
forward. 

We  now  consider  a  numerical  example.  It  is  assumed  that,  at  reconnais¬ 
sance  time  tr  =  o,  the  aggregate's  middle  point  occupied  the  position  x  *  0. 
The  fixed  numerical  assumptions  are: 

*  230  m 

*  400  m 

*  25  km/h 

*  10  km/h 

■  100  m 
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The  numerical  results  are  presented  in  figures  (8a),  (8b),  and  (8c). 

Shown  on  the  abscissa  and  on  the  ordinate  are  the  strike  time  t  (delay  between 
reconnaissance  and  strike)  and  the  probability  of  occupancy,  P  (xJ .  The 
number  of  colunns,  M,  is  odd  and  parameterized  from  1  through  7.  The  three 
graphs  refer,  in  this  order,  to  the  coordinates  =  3  km,  xt  *  4  km,  and 
x^  *  5  km. 

We  make  the  following  observations: 

1.  The  probability  of  occupancy  shows  the  typical  decline  with  increasing 
time.  This  effect  would  also  show  up  in  the  expected  number  of  tanks  killed, 
if  it  were  calculated. 

2.  The  probability  of  occupancy  increases,  in  certain  cases,  rather 
strongly,  as  the  number  of  colunns  increases.  This  effect  tends  to  soften  the 
decline  with  increasing  time. 

3.  An  increase  of  xt  has  about  the  same  effect  as  an  increase  of  t. 


53 


Figure  8. 


Probability  of  Occupation 


Figure  8a.  Aggregate  of  M  Columns,  xt  ■  3  km 


Probability  of  Occupation 


Figure  8b.  Aggregate  of  M  Columns,  *  4  km 
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Probability  of  Occupation 


Figure  8c.  Aggregate  of  M  Columns,  xt  *  5 
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0.  THE  TWO -DIMENSIONAL  TANK  ARRAY 


By  "two-dimensional  array",  we  mean  two  things.  First,  the  tank  array  has 
two  dimensions.  Second,  it  moves  in  two  dimensions.  We  employ  a  Cartesian 
coordinate  system  (x,  y)  with  both  axes  parallel  to  the  ground  and  the  x-axis 
pointing  toward  the  FEBA  (see  figure  9).  We  assume  that  the  tank  array  has  a 
rectangular  shape  and  that  its  sides  remain  at  all  times  parallel  to  the 
coordinate  axes.  At  reconnaissance  time,  t  *  o,  the  middle  point  shall 
occupy  the  position  x  *  o,  y  a  o. 

With  regard  to  the  motion  of  the  array,  we  assume  that  the  velocities  vx 
and  Vy  obey  the  independent,  Gaussian  VDFs 


(15-la) 


(15-lb) 


Only  in  x-direction  shall  the  array  have  a  systematic  velocity,  that  is 
vx  »  <  vx  >  >o  (15-2a) 

In  y-direction,  we  assume 

vy  3  <  vy  >  ■  °  (15-2b) 
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h 

\\ 


and 


x  (t)  3  < x  (t)>  3  tv 


y  (t)  3  <y  (t)>  3  tv 


(15-5a) 


(15-5b) 


For  the  probability  of  occupancy,  we  turn  to  figure  9,  where  L  and  W 
denote  the  half-length  and  the  half-width  of  the  array.  The  figure  shows  a 
rectangle  with  the  sides  parallel  to  the  coordinate  axes.  The  side  lengths  in 
x  and  y-direction  are  2L  and  2W,  respectively.  The  center,  T,  is  the 
contemplated  target  point.  It  has  the  coordinates  xt  and  yt. 

Necessary  and  sufficient  for  the  occupancy  of  T  is  that  the  middle  point 
of  the  array  lies  within  the  rectangle.  Therefore,  the  probability  of 
occupancy  is 


?o<T>  '  Gx  <xt>  Gv  <A> 


with 


x.  ♦  L 


G  (x 
x 


t  r 

t)=  j 


gx  (x)  dx 


(15-6) 


(15-7a) 


xt  ’  L 


A  Y 

Gy  (yt)  •  j  9y  (y)  dy 

yt  - L 


(15-7b) 


Finally,  let  w  denote  the  "interior  area  density"  of  tanks.  If  the  array  con¬ 
tains  N  tanks,  we  have 


N_ 

LW 


(15-8) 
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APPENDIX 


The  purpose  of  this  appendix  is  to  demonstrate  the  validity  of  the  approx¬ 
imation  implied  by  equations  (14-3)  and  (14-9).  In  order  to  save  time  and 
space,  we  prove  the  assertion  only  for  the  original  3DF,  F  (s),  rather  than 
the  RSDF,  F^  (s).  In  the  latter  case,  the  proof  would  require  the  numerical 
solution  of  some  integrals.  However,  it  is  to  be  expected  that,  If  The  approx 
imation  is  valid  for  the  SDF,  it  is  valid  also  for  the  RSDF,  though  perhaps 
not  to  the  same  degree. 

Using  now  the  SDF  instead  of  the  RSDF,  we  write  the  assertion  from  section 
N  in  the  following,  more  elaborate  form: 


+  00 


r+» 

/  ♦(Sj*  S2)  F  (ST)  F  (S2)  ds2 


dsj  ( A- 1 ) 
Sj  =  const 


•foo 

J  1  fd(2s)  Exd  . 

r~  (s  -  s0)z 

L  2  a/  j 

(A-2) 


I  S  J 


(A-3) 


For  the  proof,  we  perform  in  equation  (1)  the  substitution 


S1  +  S2 


2S ' 


S,  -  S, 


2S ' 


( A-4) 
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— 


which  is  equivalent  to 


S‘l  +  S‘2  =  S1 


S'l  “  S ' 2  S2 


(A-5) 


The  Jacobian  is 


3  (Sp  S2) 

3  (S'lf  S'2) 


=  -2 


(A-5) 


Therefore,  equation  (1)  assumes  the  form 


+oo  —op 


*  (2S1  x)  F  ( S’  x  +  S'2)  F  (S'  j  -  S ds'2 


?ds‘i 


S'  »  const 
1 


Now  observing  equation  (13-14),  we  have 


)  F  (S^  -  s*2) 


_1 _  -Q  (S^,  S'2) 

o2 s  2  *e 


(A-7) 


(A-8) 
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with 


Q  (s^,  S'2)  = 


(sysys0)2  *  (sysys/ 


2a 


2 


( A-9) 


Here,  Q  can  be  brought  into  the  following  form 


Q  (S'p  S'2)  = 


(S{  -S0)"  + 

- zz — 


( A— 10) 


If  we  heed  equation  (10)  in  equation  (7),  we  see  that  the  double-integral 
splits  into  two  single  integrals.  We  express  this  by  writing 


I 


( A— 11 ) 


where  we  still  have  some  freedom  in  defining  1^  and  I2.  We  also  nofe  that 
the  subscripts  and  the  prime  at  the  symbols  S’^  and  S’?  are  now  dispens¬ 
able.  Dropping  these  insignia,  we  write 


1 - 

o*  2 IT 


+O0 


<t>  <2S)  Exp 


(s  -  s0r 


ds 


(A-12) 


( A- 13) 
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Here,  I2  is  solvable  in  closed  form.  Using  equations  (4-4a),  (4-4b), 
and  (4-5a),  we  get 


h  '  asV*_' 

If  then  we  combine  equations  (11),  (13),  and  (14),  we  obtain 


(A-14) 


I 


1 _ 

asypP 


+CO 

r  2  'l 

1  4>  ( 2  S)  Exp  . 

(S-S0) 

- 2 — 

J 

-  oo 

°s 

s  > 

>  ds 


(A-15) 


This  equation  is  still  exact.  It  becomes  an  approximation  if  we  replace 
*s  by  as  JT.  In  this  case,  equation  (15)  goes  over  into  equation  (2), 
that  is,  I  is  approximated  by  J. 

It  is  now  in  order  to  say  a  few  words  about  the  reasons  for  and  the  accu¬ 
racy  of  the  approximation.  The  need  for  it  can  be  seen  from  equations  (14-12) 
and  (14-13).  In  equation  (14-12),  the  conditional  probabilities  of  occupation 
for  all  columns  are  summed  up.  In  equation  (14-13),  this  sum  is  then  inte¬ 
grated  with  one  and  the  same  RSDF,  FR( S) .  This,  of  course,  saves  a  lot  of 
numerical  integrations.  But  it  introduces  an  error. 

Turning  now  to  the  accuracy,  we  have  to  say  a  few  words  about  the  aimpoint 
selection.  Clearly,  the  optimal  aimpoint  is  always  the  point  where  the  proba¬ 
bility  of  occupancy  has  a  maximum.  This  means  that  the  first  bomb  should  be 
aimed  at  the  aggregate's  middle  point.  At  this  point,  the  probability  of  occu 
pancy  is  largest  for  column  3  (see  figures  7a,  7b),  somewhat  smaller  for 
columns  2  and  4,  and  still  smaller  for  columns  1  and  5.  However,  the  error 
introduced  by  the  approximation  refers  to  columns  1  and  5.  In  other  words, 
the  error  refers  to  the  smallest  contributor  to  the  probability  of  occupancy. 
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